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FAST  ALGORITHM  OF  ZONAL  WAVEFRONT  RECONSTRUCTION 
Wang  Yingbo 

Beijing  Institute  of  Technology 

ABSTRACT  This  paper  describes  u  fast  algorithm  -sub-block  wave-front  reconstruction.lt 
has  characteristics  of  smaller  required  memory  capacity,  smaller  noise  coefficient  and  better  nu 
merical  value  stabilitv.  Because  of  new  errors  of  the  algorithm,  qualitative  analysis  and  comput-Lr 
simulation  results  indicate  that  its  reconstruction  precision  is  basically  equal  to  that  of  the  usu 
algorithm. 

KEY  WORDS  adaptive  optics,  wave  — front  reconstruction,  zonal  method. 

I .  Introduction 

One  of  the  key  problems  in  adaptive  optics  is  how  to  obtain 
the  data  reconstruction  wavefront  from  the  wavefront  detector. 
Commonly  used  reconstruction  algorithms  include  the  zonal 
reconstruction  method  and  the  mode  reconstruction  method.  In  the 
zonal  reconstruction  method,  the  Southwell  approach  [1]  is 
suitable  for  a  Hartmann  detector  with  higher  reconstruction 
accuracy.  Here,  the  measured  inclination  rate  and  the  phase 
point  to  be  estimated  are  on  lattice  points .  The  phase  solution 
can  be  summarized  in  the  following  linear  algebraic  equation  sets 
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In  the  equation,  S  is  the  inclination  vector  for  two  orthogonal 
directions;  $>  is  the  phase  to  be  estimated;  D  is  a  sparse  matrix; 
and  A  is  the  matrix  related  to  the  algorithm  and  the  lattice- 
point  layout.  Eq.  (1)  is  a  constant  equation.  The  minimum  norm 
solution  of  the  least-squares  of  $  is 


A<&=DS 

(1  ) 

If  we  let 

B=A"d  ,  then 

<b  =  A~  DS 

(2) 

4>  =  BS 

(3) 

A'  is  the  Moore- 

-Penrose  inverse  of  A  in  the  broad  sense. 

After 

determining  the  system  parameters,  B  can  be  determined  off-line 
for  storage. 

If  the  number  of  the  subapertures  is  NxN,  then  B  is  the 
N2x2N2  matrix.  The  total  computational  volume  in  computing 
Eq.  (3)  is  2N4  times  the  multiplication  and  addition  operations, 
with  storage  of  2N4  elements .  Since  the  computational  volume  and 
storage  are  relatively  high,  it  is  more  difficult  to  conduct 
real-time  wavefront  reconstruction.  Based  on  the  features  of  the 
zonal  method,  the  article  presents  a  method  of  fast  calculation 
of  $ . 

II.  Subblock  Wavefront  Reconstruction  Method 

The  principle  of  the  zonal  reconstruction  method  is  to  apply 
the  measured  inclination  rate  (not  exceeding  four  points  at  the 
maximum)  of  adjacent  positions  of  the  phase  point  to  be  estimated 
for  interpolation  to  obtain  the  wavefront  phase.  If  it  is 
assumed  that  the  post-subblock  boundary  effect  of  the  subaperture 
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array  can  be  approximately  neglected,  the  wavefront 
reconstruction  carried  out  independently  for  various  blocks  with 
the  addition  of  the  determined  revised  value  is  equivalent  to  the 
above-mentioned  conventional  algorithm.  Thus,  the  computational 
volume  can  be  reduced  without  bringing  about  a  greater  effect  on 
the  estimation  accuracy.  This  is  the  fundamental  concept  of 
subblock  wavefront  reconstruction. 

Let  us  assume  that  there  is  a  square  subaperture  array  with 
NxN  subapertures,  which  are  divided  into  several  square 
subblocks.  Let  k  and  n  indicate,  respectively,  the  number  of 
subblocks  in  the  horizontal  direction,  and  the  number  of 
subapertures  in  the  horizontal  direction  in  each  block; 
obviously,  N=nk. 

After  subblocking,  the  wavefront  reconstruction  equation  of 

various  blocks  can  be  written  as: 

Aft,,-  D„s,j  ■  ^ 

4>„=AftuS‘<  =5'A 

(i  =  1  ,  2  ,  •••  .  k.  ,  j-  1,2,  •••.!')  (5) 

Since  the  sizes  of  the  subdivided  blocks  are  homogeneous  and 
equal  to  each  other,  Bij  of  the  various  blocks  are  the  same. 


Therefore,  let 


fl)  =  GS,, 


The  phase  values  obtained  from  Eqs.  (3)  and  (7)  are  not 
consistent.  This  is  because  both  solutions  to  the  equations  have 
the  minimum-norm  feature.  In  other  words,  the  mean  value  of  the 
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phases  is  zero.  After  subblocking,  this  approach  will  lead  to  a 
difference  of  a  constant  between  the  subblocked  reconstruction 
phase  and  the  reconstruction  phase  under  the  conventional 
algorithm.  Therefore,  to  ensure  continuity  of  the  wavefront,  *±j 
m  Eg.  (7)  should  have  a  revision  value  added  to  it,  which  should 
be  the  difference  between  the  reconstruction  phase  value  of  a 
specific  point  in  the  conventional  algorithm,  and  the  subblock 
reconstruction  phase  value  at  the  point.  In  other  words, 

4=0  (z\ /)-«>..  (8) 

In  the  equation,  i (k-1 ) <i ' <ik,  j (i-1 ) < j ' < jk. 

The  point  (i',  j')  should  be  located  at  the  center  position 
(n  is  an  odd  number)  or  near  the  center  position  (n  is  an  even 
number)  of  the  subblock.  Thus,  since  the  number  of  related 
points  of  the  adjacent  zones  is  unaffected  by  subblocking,  the 
phase  estimation  accuracy  at  this  point  is  not  affected. 

From  Eg.  (3),  we  know  that  $(i',  j')  in  Eq .  (8)  indicates 

the  specific  weighting  summation  of  vector  S.  In  other  words, 

4>(i'J')=bT  (i'J')S  (9) 

On  the  same  line  of  reasoning,  from  Eq.  ( 7 )  we  know  that  $±j 

( ■*"  '  ^  )  can  be  used  to  indicate  the  specific  weighting  of  vector 
Sij .  In  other  words  , 

MO) 

Since  Stj  is  a  subset  of  S,  Eq.  (10)  can  be  rewritten  as: 

mi, 

In  the  equation,  l],  (/',./')  corresponding  to  the  Si;j  portion  of  the 
S  vector  is  the  element  of  ^  ,  the  others  are  equal  to 
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zero . 


If  we  let 


hl=  b1  (f'j  ')_/t  {i’j') 


then  ..  d,rhlS 

If  the  following  equations  are  introduced, 


(12) 
(  13  ) 


~du~ 

~K- 

r  1  -■ 

dn 

H  = 

K 

/= 

i 

~dkk— 

k2x  1 

-K:  - 

k2  X  2 N2 

- 1  - 

Then  Eq.  (13)  can  be  rewritten  as 

d=  HS  f  1 4  ) 

After  revision,  the  phase  value  can  be  indicated  as 

•  I 

0=1,  2,-  ,  k  ,  1,2,-,  k)  <  ’■ 5  1 

With  respect  to  the  adaptive-optics  system  applying  the 
zonal  method,  since  the  subapertures  and  the  drivers  correspond 
to  each  other,  need  not  be  calculated.  It  is  required  only  to 
add  a  specific  DC  bias  to  all  drivers  in  the  corresponding 
subblock  after  obtaining  d^ . 

Then  the  subblock  wavefront  reconstruction  can  be  summarized 
as  requiring  the  calculation  of  Eqs .  (7)  and  (14).  If  it  is 
assumed  that  the  ratios  between  the  computational  volume  and  the 
storage,  on  the  one  hand,  and  the  corresponding  values  in  the 
conventional  algorithm  Eq.  (3)  are,  respectively,  p  and  q,  then 
it  is  easy  to  obtain  the  following: 
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1 

k~ 


1 

k: 


+  — - 
«■ 


(16) 


<7  = 


k1  ■  In'+lk'N 
2AM 

2/74+2A-:A,: 

2AM 


(17) 


For  a  given  N,  by  selecting  equal  values  or  close  to  equal 
values  between  k  and  n,  the  value  of  p  can  be  minimized.  From 
Eqs .  (16)  and  (17),  we  know  that  when  N  is  greater,  the 
computational  volume  and  the  storage  of  the  subblock 
reconstruction  method  are  reduced  by  more  than  under  the 
conventional  algorithm.  For  example,  when  N=20,  by  using  k=4  and 
n=5,  then  p  is  approximately  equal  to  0.1  and  q  is  approximately 
equal  to  0.4.  Only  when  N  is  greater  than  or  equal  to  20,  are 
the  reduced  amplitudes  one  order  of  magnitude  or  more. 


III.  Error  Analysis  and  Computer  Simulation 

With  respect  to  the  subblock  wavefront  reconstruction 
algorithm,  now  the  problem  of  greatest  concern  is  what  the 
accuracy  is.  In  the  wavefront  reconstruction  process,  there  are 
two  primary  sources  of  error.  One  is  the  sampling  error  because 
of  a  limited  number  of  subapertures,  and  the  other  is  the 
measurement  error.  The  effect  on  the  phase  estimation  accuracy 
can  be  evaluated  by  using  the  noise  propagation  coefficient, 
which  is  defined  as  the  ratio  between  the  phase  square  error  and 
the  measured  square  error  of  the  phase  difference  [1].  The 
expression  equation  is 
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is  the  element  of  matrix  B  in  Eq-.  (3).  h  is  the  spacing 
between  subapertures.  In  the  following,  the  derivation  and 
calculations  of  the  noise  propagation  coefficient  in  the  subblock 
reconstruction  method  are  carried  out.  First,  the  following 
symbols  are  introduced: 


r  vr 

; 

j~s. 

5  = 

5I: 

n  ~  x  2A':  ■ 

-  i'u- 

"CD,  ~ 

<&'= 

cD:: 

H'= 

H .. 

Nz  x  1 

-Mu- 

From  Eqs .  (13)  and  (15),  we  have 

dn  •  /=  H,/S  (19) 

Rearrange  vector  S  so  that  it  becomes  S'.  Moreover, 

corresponding  interchanges  of  various  roles  in  the  Hij  matrix  are 
carried  out  so  that  the  product  of  the  converted  H'  subject  and 
vector  S'  is  the  same.  In  other  words, 

du  ■  I  =  H'f  S '  (20) 

From  Eqs.  (7)  and  (15),  we  know  that  the  revised  phase  value 
can  be  expressed  as: 


4'  = 


r  G 


o 


5'  -H'  S‘ 


G- 

G'  S'  +  H’S' 


(21  ) 
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Let 


B  ’  =G'-H ' 
<t>'  =  B'S' 


(22) 

(23) 


Then  we  have 

Assume  that  the  noise  propagation  coefficient  in  the  subblock 
reconstruction  method  is  Ck(N),  then  from  Eq.  (18)  and  the 
foregoing  equations,  we  have 

C*(A  (Nh  (24) 

After  computation,  the  variation  status  of  Ck(N)  with  N  and 
k  is  shown  in  Fig.  1.  in  the  figure,  k=l  is  equivalent  to  the 


OS 
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0  2 
00 

Fig.  I  The  curves  of  the  noise  propagation  coefficient 

situation  in  the  conventional  algorithm.  From  the  figure  we  can 
see  that  the  propagation  coefficient  in  the  subblock 

reconstruction  method  should  be  smaller  than  in  the  conventional 
algorithm. 

The  number  of  conditions  for  solving  Eq.  (3)  can  reflect  the 
effective  degree  on  the  phase  solution  with  perturbation  in  the 
B— matrix  and  the  inclination  rate.  The  effective  degree  is 
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defined  as: 


cond  (N)  =omiJomm  (25) 

In  the  equation,  ct^  and  amin  are,  respectively,  the  maximum  and 
the  minimum  singularities  in  matrix  B.  if  it  is  assumed  that  the 
number  of  conditions  for  solving  Eq.  (7)  for  subblocks  after 
subblocking  is  the  condition  cond*  (N),  then  we  can  obtain 

cond,.  (  N  )  =  cond  (n  ) 

=  cond  (A '  k)  (26) 

Since  the  number  of  conditions  increases  rapidly  with 
increased  variation,  the  number  of  conditions  for  solving  the 
equation  is  reduced  after  subblock  reconstruction,  then  the 
stability  of  the  values  can  be  enhanced.  Fig.  2  shows  the 


Fig.  2  The  curve  of  the  condition  number 

variation  of  cond*  with  N  and  k;  (k-1)  is  equivalent  to  the 
situation  in  the  conventional  algorithm.  As  indicated  by 
analyzing  two  aspects  of  the  noise  propagation  coefficient  and 
the  number  of  conditions  mentioned  above,  it  appears  that  the 
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subblock  reconstruction  should  lead  to  smaller  wavefront 
reconstruction  errors.  However,  on  the  other  hand,  the 
consequence  of  subblocking  will  reduce  the  correlated 
reconstruction  points  of  the  phase  point  of  the  block  and  at  the 
block  boundaries.  To  some  extent,  this  will  lead  to  greater 
reconstruction  error. 

To  examine  the  effect  of  the  two  factors  mentioned  above  on 
subblock  reconstruction  accuracy,  the  author  and  his  colleagues 
conducted  computer  simulation.  During  simulation,  a  square  zone 
with  half  a  side  equal  to  1  serving  as  the  aperture,  and  the 
numbers  of  subapertures  along  the  horizontal  direction  are, 
respectively,  N=8  and  N=6;  along  the  horizontal  direction,  the 
number  of  subblocks  k=2 .  Let  us  assume  that  the  wavefront  phase 
developed  by  the  Zernike  polynomial  is  the  summation  of  the 
aperture  mode  of  26  previous  orders  [2]. 


<1>  (  .Y  -  V  )  =  X  a.  - ( -v  •  y ' 


Due  to  the  symmetry  of  aberrations  and  origins  of  various 
orders,  we  can  only  calculate  the  reconstruction  mean-square 
error  of  the  subblocks  at  the  upper  left-hand  corner;  the 
simulation  procedures  and  methods  are  shown  in  the  following 
equations . 


S>-~k--(x,y) 


s-~=J(x,y) 


( 2S ) 
(29) 
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ii 

r_ — , 

yi 

(30) 

0n 

(31  ) 

to 

H  - 

II 

-Sf 

(32) 

(33) 

aj~  V  Var  (Z  ( .v  ,  y)  -<D„) 

' - 

II 

O 

,  25 

(34) 

In  the  equations,  Sx  and  Sy  are,  respectively,  the  mean 
inclination  rates  upward  of  the  x-  and  y-directions  in  the 
subapertures .  By  using  the  mean-square  error  during 
reconstruction  in  the  conventional  algorithm,  we  can  calculate 
with  the  same  reasoning. 

The  simulation  results  are  summarized  as  follows:  . 

(1)  With  respect  to  translation,  inclination,  defocussing, 
astigmatism,  and  other  low-order  aberrations,  the  reconstruction 
mean-square  errors  in  both  algorithms  are  approximately  zero. 

This  is  because  these  wavefronts  can  be  precisely  reconstructed 
with  the  Southwell  approach. 

(2)  With  respect  to  the  coma  aberration,  and  the  vast 
majority  of  aberration  modes  higher  than  coma  aberration,  the 
reconstruction  mean  square  errors  of  both  algorithms  are 
approximately  equal  to  each  other  (with  an  accuracy  to  one- 
thousandth  of  a  wavelength) . 

(3)  In  the  aberration  modes  in  the  previous  26  orders,  there 
are  only  six  modes  with  subblock  reconstruction  mean-square  error 
relatively  increasing  over  the  case  with  the  conventional 
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algorithm.  However,  the  increasing  proportion  does  not  exceed 
5%,  or  the  absolute  values  are  smaller.  Table  1  lists  the 
simulation  results  of  these  six  modes. 


able  I  The  rms  error  of  wavefront  reconstruction  (unit  :  wavelength  ) 


^  aberration 

1 

— 

- — , 

algorithr 

n 

Z!? 

z.. 

Z; 

N  =  $ 

normal 

0.005 

0.251 

0.251 

0.033 

0.033 

0.556 

. 

in  groups 

0.025 

0.261 

1  0.261 

0.046 

0.046 

0.561 

iV  =  6 

normal 

0.008 

0.424 

0.424 

0.054 

0.054 

0.922 

in  croups 

0.043 

0.442 

0.442 

0.079 

0.079 

0.932, 

(4)  In  both  algorithms,  the  reconstruction  mean-square 
errors  increase  synchronously  with  a  reduction  in  n.  This 
indicates  that  the  magnitude  of  N  (the  wavefront  sampling 
density)  is  the  main  factor  determining  the  reconstruction  error. 

The  above-mentioned  applies  to  a  case  in  the  absence  of 
noise.  Furthermore,  simulation  computations  are  also  carried  out 
by  adding  white  noise  of  gaussian  distribution  to  the  inclination 
data.  The  signal-to-noise  ratio  (S/N)  are,  respectively,  taken 
as  30,  20,  and  lOdB.  when  S/N=20dB,  other  than  the  six  above- 
mentioned  aberration  modes  with  subblock  reconstruction  mean 
error  slightly  greater  than  that  in  the  conventional  algorithm, 
most  other  aberration  modes  indicate  a  phenomenon  with  smaller 
values  than  the  latter.  This  is  because  the  propagation 
coefficient  of  subblock  reconstruction  noise  is  smaller  and  the 
value  stabilization  is  good.  When  S/N=30dB,  the  situation  is 
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close  to  that  in  the  absence  of  noise.  The  situation  is  in  the 
mid-range  when  S/N=20dB. 

As  indicated  by  the  Zj(x,y)  linear  summation  expression 
equation  for  wavefronts  in  Eq.  (27),  and  the  above-mentioned 
computer  simulation  results,  reconstruction  accuracy  in  the 
subblock  reconstruction  algorithm  is  basically  as  that  in  the 
conventional  algorithm;  the  accuracies  may  possibly  be  slightly 
higher  in  the  presence  of  intensive  noise. 

IV.  Discussion 

As  pointed  out  by  Southwell  [1],  with  respect  to  high-order 
aberrations  in  wavefronts,  the  mathematical  model  proper 
determines  that  a  better  reconstruction  effect  cannot  be 
obtained.  However,  this  result  is  determined  by  the  value  of  N, 
the  number  of  subapertures.  Affected  by  atmospheric  turbulence 
perturbations,  the  power  spectrum  of  the  wavefront  phase  rapidly 
decreases  with  increasing  spatial  frequency  (this  also  can  be 
considered  as  the  number  of  orders  of  the  aperture  modes)  [2]. 
Therefore,  the  low-order  apertures  are  in  the  predominant 
position  in  the  wavefront  phase.  As  indicated  in  the  computer 
simulation  results,  with  a  higher  number  of  orders  of  the 
aberration  mode,  the  reconstruction  errors  in  the  subblock 
reconstruction  method  and  in  the  conventional  algorithm  increase 
synchronously  with  not  much  difference  in  the  increasing 
amplitude.  Here,  both  errors  can  only  be  reduced  by  increasing 
the  number  of  apertures.  Since  the  number  of  subapertures  is 
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limited,  the  wavefront  detector  can  detect  only  the  wavefront 
error  in  the  limited  aberration  modes.  The  foregoing  discussion 
can  explain  that  the  major  errors  are  unable  to  appear  in  the 
subblock  reconstruction  method  in  various  situations  (especially 
with  high-order  aberrations )  and  the  properties  in  the 
conventional  algorithm. 

What  has  been  discussed  in  the  preceding  sections  has  been 
the  evenly  divided  reconstruction  situation  of  N=nk.  In 
practice,  the  method  can  be  extended  to  the  situation  when  N  is 
an  arbitrary  positive  integer .  Only  by  dividing  N  into  the 
summation  of  several  numbers  close  to  each  other,  and  slightly 
revising  the  subblock  reconstruction  algorithm,  can  fast 
wave-^ront  reconstruction  still  be  carried  out,  but  only  somewhat 
increased  in  the  storage  volume  over  that  in  the  evenly  divided 
situation.  For  example,  N=17;  this,  can  be  subblocked  based  on 
N=4+4+4+5 .  Here,  two  types  of  matrix  G  in  Eq.  (7)  should  be 
stored  in  advance. 

The  Fried  method  and  the  Hudgin  method  are  also  zonal 
wavefront  reconstruction  methods  [1],  The  main  difference  from 
the  Southwell  approach  is  the  different  matrix  B  in  Eq.  (3). 
Therefore,  the  subblock  wavefront  reconstruction  method  in  this 
article  is  easily  extended  to  the  two  above-mentioned  approaches; 
close  results  are  likely  to  be  obtained. 

There  is  always  some  time  delay  in  the  adaptive-optics 
system;  this  is  one  of  the  main  factors  affecting  its  dynamic 
properties.  However,  there  is  quite  a  sizable  proportion  in  the 
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total  time  delay  when  using  the  wavefront  reconstruction.  The 
subblock  wavefront  reconstruction  can  effectively  reduce  this 
portion  of  time,  thus  enhancing  the  effect  for  wavefront 
correction,  and  the  response  rate  of  the  overall  system. 

V.  Preliminary  Conclusions 

In  the  article,  the  properties  of  the  zonal  method  are  the 
basis  for  presenting  a  fast  algorithm,  the  subblock  wavefront 
reconstruction  method.  The  problems  of  computational  volume, 
storage  quantity,  and  reconstruction  accuracy  are  analyzed  with 
computer  simulation.  Compared  to  the  conventional  algorithm,  * 
there  are  advantages  of  greatly  reducing  computational  volume, 
smaller  storage  requirement,  smaller  noise  propagation 
coefficient,  and  stabilization  of  values.  After  analysis,  the 
reconstruction  accuracy  is  basically  the  same  as  the  latter  case. 

The  author  expresses  his  gratitude  to  professor  Zhou 
Renzhong,  for  his  advice  in  the  writing  of  this  article 
The  first  draft  of  the  article  was  received  on  August  1,  1990; 
the  final  revised  draft  was  received  for  publication  on  October 
25,  1990. 


15 


MIXING  PROCESSES  DOWNSTREAM  OF  NOZZLE  BANK 
OF  A  SUPERSONIC  CW  HF/DF  LASER 

Sang  Fenting,  Gou  Chenzhou,  Zhuang  Qi,  Xie 
Xingbing,  Yuan  Qinian,  Huang  Ruiping,  and 
Zhuang  Cunhao 

Dalian  Institute  of  Chemistry  and  Physics 
Chinese  Academy  of  Sciences 


ABSTRACT  The  influence  of  the  design  parameters,  the  structure  of  the  nozzle  array  as  weli  as 
the  mode  of  H,  addition  on  the  output  power  are  discussed.  A  specific  power  of  130J;g  has  been 
obtained  with  a  mixing  distance  of  1.75mm  between  the  primary  and  secondary  nozzles. 
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I .  Foreword 

Since  the  development  of  the  first  continuous -wave  HF 
chemical  laser  by  Spencer  et  al.  in  1969,  due  to  its  advantages 
of  high  efficiency,  small  size,  no  need  for  external  power 
supply,  and  short  wavelength,  there  have  been  rapid  developments 
in  CW  HF  chemical  lasers.  In  1989,  ignition  was  successful  for 
the  HF  chemical  laser  device,  which  was  code-named  ALPHA,  as  a 

space-based  laser  weapon  in  the  United  States.  Its  power  reached 
2MW. 
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The  combustion-driven  CW  HF  (DF)  chemical  laser  device 
mainly  includes  five  components  (Fig.  1):  gas  intake  system. 


Fig.  I  Schematic  of  CW  HF  (DF  )  chemical  laser 
1.  gas  control  system  2.  combustor  3.  nozzle  4.  mirror  5.  laser  cavity  6.  exhaust  syiiem 


combustion  chamber,  nozzle  bank,  optical  cavity,  and  exhaust  gas 
system.  Controlled  by  the  gas  intake  system,  F2  and  D2  enter 
from  the  head  portion  of  the  combustion  chamber  to  be  mixed  and 
burned.  From  the  combustion,  high  temperature  splits  the  excess 
F2  into  F  atoms,  which  are  diluted  with  He  and  accelerated 
through  the  nozzle  to  enter  the  optical  cavity.  At  the  laser 
cavity  inlet,  the  atoms  mix  with  the  cavity  fuel  H2  to  form 
oscillating  excited-state  HF‘  molecules,  emitting  an  HF  laser 
beam  in  the  optical  cavity.  After  processing  in  the  gas  exhaust 
system,  exhaust  gas  is  expelled  into  the  atmosphere.  The 
properties  of  the  gas  intake  system  and  the  combustion  chamber, 
nozzle,  optical  cavity,  and  gas  exhaust  system  directly  affect 
laser  device  performance,  especially  the  nozzles.  The  properties 
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not  only  determine  the  gas  dynamic  process  of  the  entire  laser 
cavity,  but  also  has  a  direct  effect  on  laser  beam  quality. 
Therefore,  the  design,  machining,  and  experimentation  on  nozzles 
have  been  consistently  of  interest  to  researchers.  The  article 
investigates  mainly  the  effect  on  performance  of  output  power  and 
specific  power  (power  generated  per  unit  working  body  flow)  of 
the  laser  device  due  to  various  parameters  of  the  three- 
dimensional  nozzle  bank. 

II.  Nozzle  Design 

Downstream  of  the  laser  cavity  intake  (the  plane  at  the 
nozzle  exit)  there  are  the  four  following  fundamental  processes: 

(1)  flowing  out  of  the  primary  nozzles,  the  ultrasonic  gas  flow 
containing  F  and  He  atoms  has  a  diffusional  mixing  process  with 
the  H2  flow  out  of  the  secondary  nozzles. 

(2)  The  chemical  reaction  process  (generating  the  emission 

working  particles)  is  the  chemical  pump  process: 

F+  H;  —  HF'  (  v=  l  .  2 , 3  )  +  H 

(3)  excitation  and  emission  process: 

HF‘  (  r )  +  /jv  — -  HF  ‘  (  v  -  1 )  +  2/rv 

(4)  collision  relaxation  process: 

HF '  (  v )  +  A/  —  HF‘  (v-1  )  +  M 
where  M  can  be  HF,  H2,  He,  or  F2. 

In  the  foregoing  process,  the  chemical  reaction  process 
proceeds  very  rapidly  measured  in  the  order  of  microseconds . 
However,  the  diffusional  mixing  process  is  very  slow,  measured  in 
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billiseconds .  Therefore,  the  rate  of  diffusional  mixing  directly 
affects  the  amount  of  HF‘  molecules  in  the  oscillating  excited 
state.  In  other  words,  speeding  up  the  diffusional  mixing 
process  can  raise  the  laser  device  output  power. 

If  we  neglect  the  wall  thickness  at  the  nozzle  exit,  the 
diffusion  mixing  process  downstream  of  the  nozzle  exit  plane  can 
be  simplified  as  shown  in  Fig.  2. 


Fig.  2  Schematic  of  mixing  layer  structure 
1.  mirror  2.  mixing  boundary  layer 


In  the  figure,  LD  is  the  characteristic  mixing  length,  which 
is  the  axidirectional  length  when  the  cavity  fuel  H2  diffuses  to 
the  primary  nozzle  center  line.  SM  is  the  spacing  between  the 
main  and  the  secondary  nozzles  called  the  lateral— direction 


mixing  distance. 

During  the  experiments,  the  laser  cavity  pressure  pe  is 


generally  between  0.4  and  0.6kPa;  the  length  x  of  the  laser 

emission  zone  is  above  5cm.  Based  on  the  results  of 

Shackleford's  experiments  [1],  the  flowing  state  is  laminar-flow 

mixing.  The  thickness  ct(x)  of  the  laminar  mixing  level  can  be 

calculated  with  the  following  equation: 

DY 

a(x)  =  K  (  -  )'  :  ( 1  ) 

u 

In  the  equation,  K  is  a  constant  of  the  first  order  of  magnitude. 
D  is  the  diffusion  coefficient.  From  [2],  when  H2  diffuses  to 
such  gas  current  with  He  as  the  diluent,  D  is  2 . 88x1 0’A(T/p) 3/2 . 
However,  T  and  p  are  the  gas  flow  temperature  and  pressure;  u  is 
the  gas  flow  rate.  From  gas  dynamics  we  know  that  the  supersonic 

gas  flow  rate  is  w  =  7  *(7'°~r)  ;  y  is  the  specific- 

heat  ratio,  and  R  is  the  gas  constant. 

From  the  figure,  we  know  that  when  X=LD,  S(x)=SM.  Here  Eq. 

1  can  be  converted  to  LD=  K ' Si  (-^  )  .  From  the  equation,  we 

know  that  there  are  relationships  in  the  diffusional  mixing 
process  for  laminar  supersonic  flows,  the  rate  of  diffusional 
mixing  is  related  to  the  characteristic  dimensions  SM  of  the  jet 
nozzle,  as  well  as  the  gas  dynamic  parameters  p  and  T  of  the 
working  medium.  When  the  working  medium  conditions  p  and  T  are 
constant,  it  is  most  beneficial  to  accelerate  the  mixing  process 
by  reducing  the  lateral-direction  mixing  distance  SM. 

Based  on  the  above-mentioned  analysis,  the  authors  designed 
several  jet  nozzles  with  different  SM;  the  detailed  parameters 
are  shown  in  Table  1.  in  the  table,  dt  is  the  diameter  of  the 


Table  1  Elementary  parameters  of  nozzles  (calculated) 


-nozzle 

(mm ) 

I 

‘rimary  nozzle 

secondary  nozzle 

d.  (mm  ) 

n  ( No  ) 

M 

AJA, 

d t  (mm  ) 

M 

AjA, 

form  of  H;  addition 

CL-A 

4.25 

1.2 

38 

- 

25 

0.8 

1 

\ 

subsonic  fine  —  hole 

CL-n 

3.40 

1.0 

115 

5.3 

20 

0.5 

4.7 

20 

supersonic  fine -hole 

CL  - 1 

2.70 

1  0.5 

240 

:  5.3 

20 

0.6 

4.2' 

13 

supersonic  fine  — hole 

___CL  —  9 

1.75 

0.3 

450 

5.3 

20 

0.3 

3.9 

11 

supersonic  fine -hole 

jet  nozzle  throttle,  n  is  the  number  of  nozzles,  M  is  the  Mach 
number,  and  Ae/At  is  the  nozzle  exit  area  to  the  nozzle  throttle 
area  ratio.  To  find  the  lateral-direction  velocity  component  of 
the  gas  flow  and  to  increase  the  turbulence  downstream  of  the 
nozzle  exit,  two  different  nozzle  diffusion  half-angles  (11°  and 
15°)  are  considered.  Different  modes  of  hydrogen  addition  will 
have  apparent  different  gas  dynamic  phenomena  of  the  mixing 
layer.  The  authors  designed  two  modes  of  hydrogen  addition:  one 
is  hydrogen  addition  at  supersonic  speed  in  two-dimensional 


Table  2  Parameters  of  supersonic  two-dimensional  split  nozzle  H,  addition 


S* 

primary  nozzle 

secondary  nozzle 

nozzle 

(mm ) 

d,  (mm ) 

n  (No) 

M 

AJA, 

ht  (mm ) 

AJA , 

CL  -6 

1.75 

0.52 

190 

5.7 

27 

-0.2 

3.5 

CL  -7 

1.75 

0.52 

232 

5.3 

20 

0.2 

4.5 

1 

U 

1.75 

0.30 

374 

5.3 

20 

0.2 

8 

S'  =  0.0087  JT 


slits,  and  another  mode  of  hydrogen  addition  at  supersonic  speed 
in  three-dimensional  small  holes.  Table  2  shows  the  detailed 


parameters.  In  the  table,  ht  is  the  height  of  the  jet  nozzle 
when  hydrogen  was  added  to  two-dimensional  slits. 

For  convenience  in  machining,  type  H62  brass  was  used  as  the 
nozzle  material.  To  prevent  reattachment  of  the  F  atoms  onto  the 
brass  surface,  nickel  was  deposited. on  the  nozzle  inner  surface. 

Relatively  speaking,  the  boundary  layer  of  the  small- 
Reyno Id s ' -number  nozzle  is  relatively  thicker;  this  is  a  problem 
that  we  cannot  neglect.  However,  there  is  still  a  lack  of  a 
complete  effective  revision  method  on  the  additional  surface 
layer  for  nozzles  with  three-dimensional  holes.  Therefore,  the 
conventional  empirical  formula  of  the  adherent  layer  of  laminar 
flow  [3]  was  used  to  revise  the  computed  displacement  thickens  6‘ 
at  the  nozzle  exit.  The  revised  formula  is 

6'=0 . 0087x  . 

The  flow  calculation  in  the  jet  nozzle  applies  to  the 
assumption  of  one-dimensional  isentropic  flow,  regarding  r  in  the 
flow  process  in  the  nozzle  as  constant.  When  calculating  the 
mixing  parameters  and  the  exits  of  the  primary  and  secondary 
nozzles,  it  is  assumed  that  the  mixing  and  reaction  processes  are 
independent.  In  other  words,  what  first  occurs  is  the  mixing  of 
flows  from  the  primary  and  secondary  nozzles,  and  then  the 
reaction  of  the  mixed  gas  flows .  The  reaction  process  is 
adiabatic,  frictionless,  and  shock  wave-free  [4],  (This 
simplification  may  differ  somewhat  from  the  actual  flow,  but  it 
can  serve  as  a  reference. )  Tables  3  and  4  show  the  calculation 
results  for  type  CL-9  nozzles. 
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Table  3  Operating  condition  and  dimensions  for  CL  “9  nozzle 


parameter 

symbol 

units 

primary  nozzle 

secondary  nozzle 

Mach  number 

M 

5.2 

3.9 

pressure 

Pr 

torr 

5 

5 

velocity 

ut 

I  "  " 

cm  /  s 

|  2176x10-’ 

2554 x  t0: 

temperature 

r 

K 

169 

74 

density 

p 

g/  cnT 

5.51  xlO*6 

2.10  x  10’* 

mole  fraction 

yr 

0.15 

0.34 

r* 

0.51 

y  H: 

I 

displacement  thickness 

(5* 

cm 

0.0023 

0.0017 

throat  diameter 

dt 

cm 

0.030 

0.030 

exit  diameter 

dr 

cm 

0.1340 

0.1000 

expansion  half-  angle 

a 

degree 

15 

15.3 

number  of  nozzle 

n 

No 

450 

450 

The  fundamental  nozzle  structure  is  of  the  integral  type, 
without  water  cooling.  This  is  feasible  for  experiments  with 
short-term  operation.  Fig.  3  shows  the  schematic  structure. 

III.  Experimental  Results  and  Discussion 

All  thermal  experiments  with  the  nozzles  were  conducted  on 
continuous -wave  chemical  laser  device  driven  with  a  kilowatt- 
level  combustion.  In  the  experimental  process,  flows  of  the 


Table  4  Operating  condition  ( calculated  )  before  and  after  fixing  and  reaction  for  CL -9 


parameter 

before  mixing 

after  mixing 

after  reaction 

primary  nozzle 

secondary  nozzle 

M 

5.2 

3.9 

.  3.87 

1.71 

T  ( K. ) 

169 

74 

1.86 

643 

T ,  (K) 

1600 

298 

811 

1124 

349 


working  media,  F2,  D2,  He,  and  H2  were  measured;  also  measured 
were  the  combustion  chamber  pressure  pc,  laser  cavity  pressure 


Fig.  3  Schematic  of  nozzle  structure 
I.  primary  nozzle  2.  secondary  nozzle 

Pi,  distance  xc  between  the  optical  axis  and  the  nozzle  exit,  and 
the  distributions  along  the  gas  flow  direction  of  small-signal 
gain  coefficient  of  Pi=l,  2,  and  3  spectral  lines  in  the  laser 
cavity .  The  experimental  conditions  and  results  for  several  jet 
nozzles  are  shown  in  Tables  5  to  7  and  also  in  Fig.  4.  In  the 
tables,  F2/D2  is  the  molar  flow  of  F2  and  D2  for  F2/D2  entering  the 
combustion  chamber;  wx  is  the  total  working  medium  flow;  and  the 
specific  power  is  the  ratio  between  the  output  power  and  the 


total  flow  quantity  Ws. 


Table  5  Experimental  condition  and  results  for  CL  -  A  .  CL  -  11  ,  CL  -  1  and  CL -9  nozzle 


No 

s. 

(mm ) 

nozzle 

exi 

perimen 

tal  condition 

condition  of  optical  cavity 

experimental  results 

F:/D; 

(g./s  ) 

Pf 

( kg/  cm1 ) 

pL 

(torr) 

reflection 

mirror 

output 

mirror 

xr 

(mm  ) 

power 

(W) 

specific 
power  (J/g ) 

79-204 

4.25 

(CL-A) 

1.68 

8.09 

3.5 

R6M 

spherecai 

R6M 

spherecai 

20 

567 

70 

80-318 

3.40 

(CL-  11  ) 

1.63 

12.60 

1.88 

2.2 

R6M 

column 

R6M 

spherecai 

20 

1035 

82 

.  ^9-222 

2.70 

(CL- 1 ) 

1.48 

7.81 

2.30 

3.0 

R6M 

column 

R6M 

spherecai 

17 

829 

1 - 

106 

ij,  *0-359 

_  — 1 

1.75 

(CL  —  9 ) 

1.44 

7.84 

3.65 

3.0 

R6  M 

column 

R6M 

spherecai 

10  j 

1025 

130 

From  Fig.  5,  we  can  see  the  following:  (1)  when  the  lateral- 
direction  mixing  distance  of  the  primary  and  secondary  nozzles 
was  reduced  (4.55mm  to  1.75mm),  the  specific  power  of  the  laser 
device  increased  from  70J/g  to  130 J/g.  The  reason  is  due  to  the 
reduction  of  the  mixing  distance  with  accelerated  mixing  so 
that  there  is  an  increase  in  HF-  generated  per  unit  time.  This 
is  consistent  with  the  foregoing  analysis.  (2)  When  the  lateral- 
direction  mixing  distance  SM  was  reduced  and  the  mixing  was 
accelerated,  the  characteristic  mixing  length  becomes  shorter,  so 
that  the  position  of  the  optical  axis  had  to  move  forward  from 
20mm  to  10mm  from  the  nozzle  exit  so  that  the  equivalent  light 
release  zone  became  shorter.  In  Fig.  4,  as  indicated  with  the 
gain  distribution  of  the  p  branch,  there  was  a  positive  gain 


Table  6  Experimental  results  for  CL  -  9  nozzle 


No 

ex 

perirnent 

al  condition 

condition 

of  optical  c; 

avity 

experimental  resu|ls 

F;/D; 

h; 

(g/  s ) 

reflection 

mirror 

output 

mirror 

xc 

(mm ) 

power 

(W) 

specific 
Power  (j/g) 

80-354 

1.41 

7.96 

3.55  • 

3.0 

R6M 

column 

R6M 

spherical 

10 

959 

120 

80-  355 

7.84 

3.60 

3.0 

R6M 

column 

R6“ 

spherical 

10 

917 

117 

80-356 

1.41 

7.98 

3.65 

3.0 

R6M 

column 

R6M 

spherical 

10 

934 

117 

80-  357 

1.43 

7.84 

3.65 

3.0 

R6M 

column 

R6M 

spherical 

10 

933 

119 

80-  358 

:  1.44 

7.83 

3.65 

3.0 

R6M 

column 

R6M 

spherical 

10 

955 

122 

80  —  359 

1.44 

!  7.84 

3.65 

3.0 

R6M 

column 

R6M 

spherical 

10 

1025 

131 

i 

80-360  1 
1 

1.44 

7.84 

3.65 

3.0 

10 

933 

119 

coefficient  in  the  type  CL-9  nozzle  only  in  the  range  of 
approximately  2cm.  At  further  distances,  it  became  an  absorption 
region . 

From  Table  7  we  can  see  that,  with  respect  to  the  same 
nozzle  throttle  dimensions,  although  the  distance  SM  was  reduced 
for  mixing  in  the  lateral-direction,  yet  since  the  two- 
dimensional  hydrogen  addition  was  applied,  the  specific  power  was 
reduced  by  nearly  20%.  However,  for  type  CL-05  and  type  CL-9 
nozzles  with  dt=0.3mm,  the  lateral-direction  mixing  distance  was 
1.75mm,  but  the  specific  power  was  38%  less.  This  occurred 
mainly  when  there  was  hydrogen  addition  at  the  two-dimensional 


Table  7  Comparison  between  two-dimensional  split  nozzle  and  three  dimensional  fine -hole  H,  addition 


No 

nozzle 

d, 

(mm) 

i 

experimental  condition 

condition  of  optical  cavity 

experimental  results 

F:/D. 

w; 

(Sis) 

Pf 

(kg/ cm;  ) 

Pt 

(torr 

reflection 

mirror 

output 

mirror 

,T( 

(mm ) 

power 

(W  ) 

specific 
power  (J  g ) 

80-008 

1 

— j 

O 

0.5 

1.44 

8.71 

3.20 

2.5 

R6M 

spherical 

R10M 

spherical 

13 

881 

101 

80-119 

CL  —  6 

0.5 

1.42 

9.05 

3.20 

/ 

R6M 

spherical 

R10M 

spherical 

10 

721 

80 

80-097 

CL  —  7 

0.5 

1.47 

9.26 

2.65 

© 

C'\ 

R6M 

spherical 

R 1 0M 

spherical 

1  13 

730 

79 

j 

SO  -  156  j 

l 

CL -5 

0.3 

1.31 

7.85 

4.60 

2.0 

R6M 

spherical 

R1 0“ 

spherical 

10 

638 

SI 

1  1  1 
80-359  |  CL  — j  0.3 

1.44 

7.84 

3.65  . 

3.0 

R6M 

spherical 

R6M 

spherical 

10 

1025 

130 

A  Supersonic  three-dimensional  fine  hole  H.  addition. 


slits,  the  length  of  the  mixing  surface  of  the  gas  flow  from  the 
primary  and  secondary  nozzles  was  smaller  than  in  the  three- 
dimensional  case.  Moreover,  the  nozzle  expansion  ratio  was 
smaller  for  hydrogen  addition  at  the  two-dimensional  slits. 

During  the  experiment,  the  pressures  at  the  exits  of  the  primary 
and  secondary  nozzles  were  seriously  mismatched  so  that  the  gas 
dynamic  conditions  of  the  flow  field  deteriorated. 

IV.  Conclusions 

When  using  a  three-dimensional,  small-hole  nozzle  array  for 
a  combustion-driven  HF  chemical  laser  device,  it  was  beneficial 
to  upgrade  the  specific  power  of  the  HF  laser  device  to  reduce 
the  lateral-direction  mixing  distance  between  the  primary  and  the 
secondary  nozzles  and  to  adopt  hydrogen  addition  with  three- 
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dimensional  supersonic  speed.  As  indicated  in  the  experiments, 
when  the  array  SM  was  reduced  from  4.25mm  to  1.75mm,  the 
corresponding  throttle  diameter  of  the  primary  and  secondary 
nozzles  was  reduced  from  1.2mm  to  0.3mm,  and  the  specific  power 
of  the  HF  laser  device  was  upgraded  from  70J/g  to  130J/g. 

The  authors  are  grateful  to  Fang  Jingke,  Li  Mingsheng,  Min 
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